泾川流域土地利用/覆被时间序列数据集(1976‒2015)

王 荣1,3,王亚琴2*,胡良柏1

1. 天水师范学院资源与环境工程学院,天水 741001;
2. 中国国土勘测规划院,北京 100035;
3. 兰州交通大学测绘与地理信息学院,兰州 730070

摘 要:泾川流域指渭河的最大支流——泾河上游、泾川县城以上的河流汇水区。泾川流域土地利用/覆被分类体系由6个一级类与14个二级类组成,一级类包括林地、草地、湿地、耕地及人工表面等6类。基于1976–2015年多时相Landsat MSS/TM/OLI系列影像及eCognition、ENVI平台,采用决策树分类、面向对象基于知识规则集信息提取等方法,得到泾川流域时间序列土地利用/覆被变化数据集(1976‒2015)。该数据集包括:(1)泾川流域边界数据;(2)1976‒2015年10期(包括:1976、1991、1995、2000、2003、2005、2007、2010、2013、2015)土地利用/覆被分布数据,其中,1976年数据分辨率为60 m;2013、2015年数据分辨率为15 m;其他年份数据分辨率为30 m;(3)1976‒1991、1991‒2000及2000‒2015年3期土地利用/覆被变化数据,其中2000‒2015年数据分辨率为15 m;其他2期数据分辨率为30 m。数据集存储为.shp和.tif格式,由73个文件,数据量为252 MB(压缩为2个文件,数据量为6.83 MB)。

关键词:泾川流域;土地利用/覆被;1976–2015;黄河流域

DOI: https://doi.org/10.3974/geodp.2022.01.14

CSTR: https://cstr.escience.org.cn/CSTR:20146.14.2022.01.14

数据可用性声明:本文关联实体数据集已在《全球变化数据仓储电子杂志(中英文)》出版,可获取:
https://doi.org/10.3974/geodb.2021.09.03.V1https://cstr.escience.org.cn/CSTR:20146.11.2021.09.03.V1.

1 前言

泾河是黄河第一支流渭河的第一大支流,发源于六盘山东麓,位于西北黄土高原106°14′E–108°42′E、34°46′N–37°19′N,流经陕甘宁三省,注入渭河。泾河流域灌溉条件良好,农业发达,经济开发潜力大,是西北地区产粮区[1],但流域内地形支离破碎、植被覆盖度低,林草地占流域总面积10%左右,水土流失严重,因此,2000年被列为全国首批退耕还林(草)工程示范区2014年被列为新一轮退耕还林脱贫攻坚工程[2]。泾川流域指泾河上游流域泾川水文站以上的河流汇水区,位于泾河流域西南部,地势西高东低,属于半湿润过渡带,具有黄土高原暖温带气候和山地气候特点[3]。根据流域内平凉气象站多年观测资料,泾川流域多年年均气温8.9 ℃,年降水量355–845 mm,夏季降水量偏多,降水年际间变化差异显著。近年,国家开展大规模黄河流域生态恢复工程建设,退耕还林(草)、限制乱采乱伐乱开发等土地利用方式的改变,使得泾河流域土地利用/覆被发生巨大变化[4,5],从而导致泾河流域生态系统结构发生相应变化。

因此,研究结合泾川流域区域特色,借鉴生态十年环境遥感土地利用/覆被分类体系[6],以美国陆地卫星Landsat TM/OLI系列遥感影像为基础数据源,采用面向对象基于规则集分类、决策树分类等信息提取与目视解译相结合的方法,提取泾川流域1976–2015年长时间序列土地利用/覆被分布及其变化数据集,并通过Google高分辨率影像、实地采样点数据及已有生态十年(2010年)成果数据对提取的土地利用数据进行精度验证,以全面了解泾川流域土地利用/覆被类型及其动态变化趋势。

2 数据集元数据简介

《泾川流域时间序列土地利用/覆被变化数据集(1976–2015)》[7]的元数据简介见表1

表1 《泾川流域时间序列土地利用/覆被变化数据集(1976‒2015)》元数据信息表
条目描述
数据集名称泾川流域时间序列土地利用/覆被变化数据集(1976-2015)
数据集短名LanduseJingChuanBasin_1976-2015
作者信息王荣,天水师范学院,0119061@stu.lzjtu.edu.cn
王亚琴,中国国土勘测规划院,freefly_99@126.com
胡良柏,天水师范学院,398361732@qq.com
地理区域泾川流域
空间分辨率1976年,60 m×60m;1991、1995、2000、2003、2005、2007、2010年,30 m×30m;2013、2015年,15 m×15m
数据年代1976、1991、1995、2000、2003、2005、2007、2010、2013、2015
数据格式泾川流域边界数据,.shp;泾川流域土地利用/覆被类型及其动态变化数据,.tif
数据量252 MB 时间分辨率 年
数据集组成1∶50000泾河流域边界.shp格式数据;1976–2015年10期.tif格式土地利用/覆被分布数据;1976–1991、1991–2000及2000–2015年3期.tif格式土地利用/覆被类型变化数据
基金项目天水师范学院(CXJ2020-14);甘肃省教育厅(2019B-134)
数据制作工具eCognition、ArcGIS10.2、ENVI5.3
出版及共享服务平台全球变化科学研究数据出版系统 http://www.geodoi.ac.cn
地址北京市朝阳区大屯路甲11号100101,中国科学院地理科学与资源研究所
数据共享政策(1)"数据"以最便利的方式通过互联网系统免费向全社会开放,用户免费浏览、免费下载;(2)最终用户使用"数据"需要按照引用格式在参考文献或适当的位置标注数据来源;(3)增值服务用户或以任何形式散发和传播(包括通过计算机服务器)"数据"的用户需要与《全球变化数据学报(中英文)》编辑部签署书面协议,获得许可;(4)摘取"数据"中的部分记录创作新数据的作者需要遵循10%引用原则,即从本数据集中摘取的数据记录少于新数据集总记录量的10%,同时需要对摘取的数据记录标注数据来源[8]
数据和论文检索系统DOI, CSTR, Crossref, DCI, CSCD, CNKI, SciEngine, WDS/ISC, GEOSS

3 数据研发方法

3.1 数据来源

本数据集制作所采用Landsat MSS/TM/OLI系列影像数据来源于USGS(https://earthexplorer.usgs.gov/),根据土地利用/覆被解译需求,结合研究区影像云量及质量,获取1976、1991、1995、2000、2003、2005、2007、2010、2013、2015共10期21幅影像,传感器、数据时间、空间分辨率及行列号如表2;30 m空间分辨率生态十年(2010年)数据来源于中国生态系统调查与评价;辅助野外调查的30 m DEM数据来源于美国地质调查局;泾川流域平凉站逐月气象数据来源于国家气象科学数据中心中国地面气候资料月值数据集。结合信息提取结果地类面积,随机布设采样点,2013年开展了泾川流域的实地考察,获得部分验证采样点土地利用/覆被信息。部分地区受交通条件限制,主要通过Google Earth获取验证点,最终,共获得地面验证采样点177个,实地调查采样点89个,Google Earth影像布点88个,如图1所示。

表2 《泾川流域时间序列土地利用/覆被变化数据集(1976–2015)》产品数据源表
年份传感器数据时间空间分辨率轨道号
1976MSS1976-6-1060 m138035/138036
1991TM1991-8-23, 1991-8-3030 m128035/129035
1995TM1995-8-18, 1995-8-2530 m128035/129035
2000TM2000-4-16, 2000-7-3030 m128035/129035
2003TM2003-6-6, 2003-6-530 m128035/129035
2005TM2003-8-15, 2005-7-12, 2005-10-730 m128035/129035
2007TM2007-5-6, 2007-8-3, 2007-9-2030 m128035/129035
2010TM2010-4-28, 2010-5-2330 m128035/129035
2013OLI2013-10-6, 2013-11-1415 m128035/129035
2015OLI2015-5-12, 2015-7-2415 m128035/129035
实地调查采样及样点空间分布

图1 实地调查采样及样点空间分布

下载

3.2 研发流程

研发结合泾川流域区域特色,以生态十年环境遥感监测土地覆盖分类体系[6]为参考,制定泾川流域土地利用/覆被分类体系,删除林地对应的二级类中常绿阔叶林、常绿阔叶灌木林、常绿针叶灌木林及乔木、灌木园地等;泾川流域无湿地及湖泊,因此,湿地对应二级分类体系中仅保留水库/坑塘及河流;泾川流域土地利用/覆被分类体系及代码如表3。其次,基于Landsat MSS/TM/OLI系列影像,在决策树分类及基于eCognition软件平台面向对象多尺度分割的基础上,结合DEM及其衍生坡度等数据,确定各地物目标类的解译标志库,建立规则库,提取泾川流域土地利用/覆被类型,在此基础上,通过目视解译,对面向对象信息提取结果进一步处理;最后,基于Google Earth高分辨率遥感影像及实际调查采样点、生态十年数据(2010年)数据对提取的泾川流域土地利用/覆被信息数据精度进行检验与评价,以分析1976‒2015年泾川流域土地利用/覆被类型及其动态变化趋势,研发技术路线如图2所示。

泾川流域土地利用/覆被及其动态变化数据集研发技术路线图

图2 泾川流域土地利用/覆被及其动态变化数据集研发技术路线图

下载

4 数据结果与验证

4.1 数据集组成

《泾河流域土地利用/覆盖时间序列数据集(1976‒2015)》包括73个数据文件,分别为土地利用/覆被分布数据集、土地利用/覆被变化数据集、泾川流域边界数据集。由下列数据组成:

(1)土地利用/覆被分布数据集:包括1976、1991、1995、2000、2003、2005、2007、2010、2013、2015年共10期土地利用/覆被类型分布数据(图3),.tif格式;时间分辨率为年,空间分辨率1976年为60 m×60m,1991、1995、2000、2003、2005、2007、2010年为30 m×30m,2013、2015年为15 m×15m;以"LandCover_JingchuanBasin_年份.tif"形式命名,其中栅格值代表土地利用/覆被类型代码,具体见表3,如栅格值"41"表示"旱地";

(2)土地利用/覆被变化数据集:包括1976–1991、1991–2000及2000–2015年共3期土地利用/覆被类型变化数据(图4),.tif格式;时间分辨率分别为16、10、16年,空间分辨率1976–1991、1991–2000年为30 m×30m,2013–2015年为15 m×15m;以"LandCoverChange_JingchuanBasin_起始年份_终止年份.tif"形式命名,如1976–1991年土地利用/覆被变化数据命名为LandCoverChange_JingchuanBasin_1976_1991.tif,其中栅格值代表土地利用/覆被类型转化,如"4121"表示土地利用类型由"旱地(41)"转化为"草地(21)";

(3)泾川流域边界数据集:".shp"格式的1∶50000泾川流域边界数据集以"Jingchuan Basin_Boundary.shp"形式命名。

1976-2015年泾川流域10期土地利用/覆被分类数据空间分布

图3 1976–2015年泾川流域10期土地利用/覆被分类数据空间分布

下载
表3 数据集使用的分类体系
Ⅰ级分类及代码Ⅱ/Ⅲ级分类代码及分类体系Ⅰ级分类及代码Ⅱ/Ⅲ级分类代码及分类体系
1 林地101 落叶阔叶林3 湿地32 河流
102 常绿针叶林4 耕地41 旱地
103 落叶针叶林5 人工表面51 居住地
104 针阔混交林52 工业用地
105 落叶阔叶灌木林53 交通用地
2 草地21 草地54 采矿场
3 湿地31 水库/坑塘6 其他61 裸土
1976-2015年泾川流域3期土地利用/覆被动态变化可视化图

图4 1976–2015年泾川流域3期土地利用/覆被动态变化可视化图

下载

4.2 数据结果

4.2.1 1976–2015年泾川流域土地利用/覆盖

通过1976–2015年10期土地利用/覆被土地利用空间分布图(图3)对比分析可知:1976年泾川流域土地利用类型相对单一,1991–2015年土地利用类型趋于稳定,1991‒2015年草地及居住地面积显著增大,针阔混交林、落叶阔叶灌木林及落叶阔叶林主要分布在高海拔地区,居住地主要分布在泾河、颉河、后峡河源头及泾河两侧的河谷地带;2003年以来采矿用地面积逐渐扩大,主要分布在平凉市崆峒区南部的岩石采矿区,而工业用地主要分布在崆峒区中部,泾河南侧的工业园区,水库主要以泾源县卧龙山水库为主。

结合表4可知:1976‒2015年泾川流域土地利用/覆被类型中旱地面积逐渐减少,减少了647.10 km2,草地、落叶阔叶灌木林面积逐年增加,分别增加了245.44、252.19 km2;人类活动用地面积,除居住地面积1991年扩张程度最大,年均扩张面积达2.12 km2/a,1991‒2015年泾川流域交通用地、工业用地、采矿场及水库/坑塘等人类活动用地面积逐年增加,且交通用地面积增加最多,2015年达17.13 km2,是1976年的18.62倍。

表4 1976‒2015年泾川流域土地利用/覆被面积统计(单位:km2
土地利用类型年份(年)
1976199119952000200320052007201020132015
X112.5142.1950.8571.8271.6979.6982.2087.9494.5895.32
X215.0317.2918.1118.1922.6822.6222.6122.9322.8022.83
X31,776.951,929.141,892.251,594.631,565.871,522.031,345.471,264.611,187.541,124.85
X443.2032.1227.2029.3220.7824.5520.8320.6718.4921.28
X5653.94469.14498.37731.18703.47725.50889.07919.10917.89899.38
X6153.76152.37162.16164.61161.77161.37160.73163.70166.30166.59
X7309.30293.83291.76328.07387.07395.34406.49442.96510.82562.50
X81.800.2119.62
X9180.46194.01186.67188.25191.34192.08191.77194.27195.80195.77
X100.929.2210.8410.6211.3012.3814.4716.2316.8917.13
X110.890.891.050.840.991.811.171.171.71
X124.1312.4113.7013.6613.7413.7513.7613.8813.8613.83
X131.131.301.662.052.842.426.89
X140.010.050.060.840.711.422.403.614.56

注:居住地(X1)、常绿针叶林(X2)、旱地(X3)、河流(X4)、草地(X5)、落叶阔叶林(X6)、落叶阔叶灌木林(X7)、裸土(X8)、针阔混交林(X9)、交通用地(X10)、水库/坑塘(X11)、落叶针叶林(X12)、工业用地(X13)、采矿场(X14);表中"–"表示两类土地利用类型间无转换。

4.2.2 1976‒2015年泾川流域土地利用/覆盖动态变化

1976‒1991年、1991‒2000年及2000‒2015年三期土地利用动态变化数据(图4)及其转移矩阵(表5-表7)。

表5 1976‒1991年泾川流域土地利用/覆被动态转移矩阵(单位:km2

转移前类型转移后的土地利用/覆被类型求和
X1X2X3X4X5X6X7X9X10X11X12
X11.810.031.84
X20.04--0.04
X319.052.321.411.405.424.8915.44
X40.1010.730.080.2511.06
X50.160.88146.782.8911.865.200.243.22170.19
X60.120.816.880.197.88
X719.822.640.850.282.7627.59
X81.190.151.34
X90.070.310.040.060.48
求和19.311.00178.522.474.125.8418.1312.935.270.536.17233.98

注:居住地(X1)、常绿针叶林(X2)、旱地(X3)、河流(X4)、草地(X5)、落叶阔叶林(X6)、落叶阔叶灌木林(X7)、裸土(X8)、针阔混交林(X9)、交通用地(X10)、水库/坑塘(X11)、落叶针叶林(X12);表中"–"表示两类土地利用类型间无转换。

通过图4a及表5可知:1976–1991年草地转化为其他土地利用类型的面积最大,总面积达170.19 km2,年均转换率达11.35 km2/a;其他土地利用类型,包括草地、落叶阔叶灌木林及河流等,转换为旱地的面积最大为178.52 km2,其中草地转化为旱地的面积达到最大,为146.78 km2;综上可知:1976–1991年泾川流域大量开垦草地,增加旱地种植面积。该阶段人类活动用地以居住地、工业用地及交通用地为主,其他土地利用类型转化为人类活动用地总面积为24.58 km2,且以转化为居住地和交通用地为主,人类活动用地转化为其他土地利用类型面积仅为3.20 km2

结合图4b与表6可知:1991‒2000年旱地转化为其他土地利用类型的面积最大,达363.93 km2,主要分布在中低海拔地区,其中旱地转化为草地的面积最大,为301.23 km2;其次为草地转化为旱地、旱地转化为落叶阔叶灌木林,面积分别为37.99、32.13 km2,即1991‒2000年泾川流域草地开垦为耕地的同时,大面积开始实施退耕还草还林。1991‒2000年以来,人类活动用地类型多样化,包括居住地、交通、采矿场及工业用地等,其他土地利用类型转化为人类活动用地的面积为27.44 km2,且以转化为居住地和交通用地为主,相反,人类活动用地转化为其他土地利用/覆被类型的面积仅为1.43 km2,即与1976‒1991年相比,人类活动用地面积净增加了26.01 km2

图4c及表7表明:与1991‒2000年的变化趋势相同,2000‒2015年旱地转化为其他土地利用/覆被类型的面积最大,达509.31 km2,其中旱地转化为草地、落叶阔叶灌木林的面积最大,分别为345.65、104.61 km2;其次为草地转化为落叶阔叶灌木林,转化面积达129.99 km2,即2000‒2015年仍实施退耕还林(草)工程,且草地转换为林地的面积有所增大,1991年以来,国家开始注重泾川流域生态恢复与保护[2]。该阶段其他土地利用/覆被转化为人类活动用地面积为64.74 km2,但仍以转化为居住地为主,人类活动面积净增加了53.8 km2,比1991‒2000年增加了27.79 km2,包括居住地、裸地、水库/坑塘、采矿用地、交通、工业用地等。

表6 1991‒2000年泾川流域土地利用/覆被动态转移矩阵(单位:km2

转移前类型转移后的土地利用/覆被类型求和
X1X2X3X4X5X6X7X9X10X11X12X13
X11.081.08
X324.102.85301.230.4732.130.711.321.12363.93
X40.317.047.35
X50.3537.990.0914.932.7656.12
X60.040.301.312.263.91
X74.5513.780.262.7621.35
X80.050.150.20
X90.741.8610.890.821.3315.64
X100.050.05
X110.100.10
求和24.850.7450.812.94317.1711.6249.198.491.320.151.331.12469.73

注:居住地(X1)、常绿针叶林(X2)、旱地(X3)、河流(X4)、草地(X5)、落叶阔叶林(X6)、落叶阔叶灌木林(X7)、裸土(X8)、针阔混交林(X9)、交通用地(X10)、水库/坑塘(X11)、落叶针叶林(X12);表中"–"表示两类土地利用类型间无转换。

表7 2000‒2015年泾川流域土地利用/覆被动态转移矩阵(单位:km2

转移前类型转移后的土地利用/覆被类型求和
X1X2X3X4X5X6X7X8X9X10X11X12X13X14
X18.250.351.230.240.070.6910.83
X330.460.192.25345.651.22104.6114.390.104.410.094.561.38509.31
X40.427.020.190.568.19
X50.590.4846.040.123.19129.993.525.110.360.260.06190.51
X61.320.240.160.280.860.133.03
X70.251.285.450.3914.220.170.292.231.5825.86
X90.030.460.240.030.76
X110.080.08
X130.030.03
求和31.723.2767.143.27362.034.58235.9418.207.575.400.380.065.253.79748.60

注:居住地(X1)、常绿针叶林(X2)、旱地(X3)、河流(X4)、草地(X5)、落叶阔叶林(X6)、落叶阔叶灌木林(X7)、裸土(X8)、针阔混交林(X9)、交通用地(X10)、水库/坑塘(X11)、落叶针叶林(X12)、工业用地(X13)、采矿场(X14);表中"–"表示两类土地利用类型间无转换。

4.2.3 1960‒2015年泾川流域气象要素逐月变化趋势

泾川流域日照时数、月最低气温、月最高气温、月平均气温如图5。图5a表明:月最低气温、最高气温及月平均气温变化趋势完全一致,呈拱型分布,1–7月气温逐渐升高,7月气温达到最大,7–12月气温开始逐渐降低,且拱桥两端月平均气温与月最低气温差较小,拱型中部6、7月的月最高温、月最低温与月平均气温的气温差最大;二次函数是拟合气温,包括月最低气温、月最高气温及月平均气温,与月份的最佳拟合函数关系,R2≥0.910,3。通过图5b可知:整体来看,1960‒2015年泾川流域日照时数平均日照时数在200h左右,月份上呈拱型分布,7月日照时数较大,2月、11月日照时数相对较少,其中2000年7月日照时数最大,达304.3 h;1975年9月日照时数最小,仅为60.2 h。1月、3–5月、12月年际日照时数分布相对集中,而2月、9–11月日照时数年际分布较分散。

1960-2015年泾川流域气象要素变化趋势统计

图5 1960–2015年泾川流域气象要素变化趋势统计

下载

4.3 数据结果验证

因数据集时间序列较长,对1976‒2015年10期土地利用/覆被信息提取结果采用地面实地调查或高分辨率影像精度检验难度大,现实可行性差,因此,研究仅对2010年及2013年数据进行了精度检验,以代表所有时间序列影像解译精度。利用2013年实地调查采样点及高分辨率遥感影像采样点对2013年影像解译结果进行精度验证,基于已有生态十年(2010年)数据及高分辨率遥感影像采样点数据对泾川流域2010年土地利用/覆被解译结果精度进行验证。从混淆矩阵看,Kappa系数≥0.8532,二级分类精度≥86.65%,一级分类精度达90%以上。一般认为Kappa系数大于0.61,分类效果较好,即研究获取的土地利用覆盖数据分类精度可靠。各土地利用类型中,水库/坑塘、居住地、裸地及工业用地等人类活动用地、草地、旱地、落叶阔叶灌木林的识别精度较高,达82%–96%之间,但针阔混交林与落叶阔叶林等有时较难区分,特别是光谱波段较接近的植被类型。

5 讨论与总结

泾河是黄河第一支流渭河的第一大支流,泾川流域位于泾川水文站以上的泾河上游流域,是退耕还林(草)工程及黄河流域生态保护与治理示范区,明确泾川流域土地利用/覆被类型及其变化趋势对泾河流域生态保护、流域规划及制定引导政策等具有重要意义。研发基于1976、1991、1995、2000、2003、2005、2007、2010、2013、2015年10期Landsat MSS/TM/OLI系列不同分辨率影像数据集,采用遥感影像分割、基于知识规则集的面向对象与传统决策树分类等信息提取方法相结合,提取了1976–2015年泾川流域土地利用/覆被类型及其动态变化信息,二级分类精度≥86.65%,并分析了1960–2015年气象要素的变化趋势,为黄河流域典型示范区国家退耕还林(草)工程及生态恢复与保护提供支撑,主要结论如下:

(1)1976‒2015年泾川流域土地利用/覆被类型中,旱地面积不断减小,草地及落叶阔叶灌木林面积逐年增加,旱地减少面积与草地及落叶阔叶灌木林地增加的面积基本一致;人类活动用地面积逐年略有增加,其中居住地1991年增加的程度最显著。

(2)尽管1976‒2015年人类活动用地面积不断略有增加,但1991年以来泾川流域土地利用/覆盖的类型基本保持稳定;1976‒1991年为大面积开垦草地为旱地的阶段,1991–2000年以退耕还草为主,2000‒2015年落叶阔叶灌木林的面积不断增加,即1991年以来国家通过退耕还林(草)、退草还林等工程,逐渐开始重视及加强泾河流域生态恢复及保护,为黄河流域生态保护和高质量发展奠定了基础。

(3)泾川流域日照时数、月最低气温、最高气温及月平均气温变化趋势基本一致,呈拱型分布,7月份均达到最大,且二次函数是拟合1960‒2015年泾川流域气温与月份间变化趋势的最佳函数,R2≥0.910,3。

作者分工:王亚琴对数据集的研发做了总体指导;王荣分析并撰写了数据论文;胡良柏收集和处理了数据。

利益冲突声明:本研究不存在研究者以及与公开研究成果有关的利益冲突。

参考文献